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We report on hyperfine-resolved spectroscopic measurements of the electric-dipole forbidden 
5P3/2 — + 8pi/2 transition in a sample of ultracold 87 Rb atoms. The hyperfine selection rules en- 
able the weak magnetic-dipole (Ml) contribution to the transition strength to be distinguished 
from the much stronger electric-quadrupole (E2) contribution. An upper limit on the Ml transition 
strength is determined that is about 50 times smaller than an earlier experimental determination. 
We also calculate the expected value of the Ml matrix element and find that it is less than the 
upper limit extracted from the experiment. 

PACS numbers: 32.70.Cs, 31.10.+Z, 32.80.Ee 



I. INTRODUCTION 

The dominant electronic transitions in atoms are gen- 
erally electric dipole (El) transitions, which for single- 
electron atoms obey the selection rule A£ = ±1, where I 
is the electron orbital angular momentum. The resonance 
lines in alkali atoms (e.g., 5s— >bp in Rb), with oscilla- 
tor strengths on the order of unity, are prototypical ex- 
amples. Obviously, there are many potential transitions 
that are not allowed by the El mechanism. The lead- 
ing order electric-dipole-forbidden transitions are electric 
quadrupole (E2) and magnetic dipole (Ml) transitions. 

Although these dipole-forbidden transitions are gener- 
ally rather weak, they can in fact be quite useful and 
important. On the practical side, they expand the num- 
ber of states that can be accessed by single-photon or 
multi-photon laser excitation. We have recently exploited 
this to probe high-lying d states of Rb by single-photon 
E2 transitions from the 5s ground state [1[. Forbidden 
transitions can provide narrow resonances and thus play 
an important role in frequency standards and atomic 
clocks. In this context, optical-lattice-induced frequency 
shifts due to such transitions must also be considered [2| • 
States that can decay only by dipole-forbidden transi- 
tions are usually metastable and thus play an important 
role in accumulating population in excited gases such as 
discharges. On the more fundamental side, these El- 
forbidden transitions can be sensitive probes of atomic 
wavefunctions and thus serve an important role in testing 
atomic structure calculations. Such tests are important 
to the interpretation of atomic measurements of funda- 
mental symmetries such as parity violation [1, Q. Ml 
transitions between states of different principal quantum 
number n are especially useful in this regard as they are 
nominally subject to the An=0 selection rule. Viola- 
tions of this rule are due to relativistic and spin-orbit 
effects. Measuring these Ml transition strengths in var- 
ious atomic systems and comparing them to theory can 
help determine the relative importance of contributions 



such as negative-energy states @. 

In the present work, we explore the El-forbidden bp — * 
8p transition in Rb by laser excitation of ultracold atoms. 
The use of ultracold atoms virtually eliminates Doppler 
broadening and allows hyperfine-selective state prepara- 
tion Q. With this ability to resolve specific hyperfine 
transitions, we exploit the associated hyperfine selection 
rules to separate the E2 and Ml contributions to this El- 
forbidden transition. Specifically, we use 87 Rb, with a nu- 
clear spin 1=3/2, and prepare atoms in the 5p3/2(-F /= 0) 
state. As shown in Fig. 1, and derived in Sec. UH the 
transition to Spi/ 2 (F"=1) is allowed only by the Ml pro- 
cess, while the transition to 8p 1 / 2 (i 7 "'=2) is only E2 al- 
lowed. 

Our work was motivated in part by the surprising re- 
sults of Bayram, et al. p} . They used a polarization tech- 
nique to measure the relative E2 and Ml contributions to 
the Rb bp — > 8p transition and obtained an anomalously 
large Ml fraction, at odds with theoretical calculations. 
Here we resolve this discrepancy in favor of theory, set- 
ting an upper limit on the Ml contribution which is larger 
than, and therefore consistent with, the predicted value. 

There have been a number of other E2/M1 ratio mea- 
surements in various atomic systems. For example, in 
atomic bismuth, fits to the hyperfine spectra of light 
emitted by a discharge revealed the relative importance 
of E2 and Ml transitions @,|^. The M1-E2 interference 
effect in the Zeeman spectra has been utilized in atomic 
lead [13, [HI- I n atomic thallium, measurements of ab- 
sorption spectra and Faraday rotation, combined with 
lineshape modeling, yielded the E2/M1 ratio 12]. In all 
of these examples, the Ml process was dominant because 
the transitions took place within the same ground-state 
electronic configuration. In our case, both initial and 
final states are electronically excited and the transition 
changes n, so the Ml process is significantly suppressed. 
Our work also differs in that the E2-only and Ml-only 
transitions are well resolved in our excitation spectra, 
obviating the need for sophisticated modeling. We note 
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that there have also been proposals for measuring E2/M1 
ratios in alkali n\p — n 2 p transitions using polarization 
effects in two- wave mixing [l3|, ■ 

This paper is organized as follows. In Sec. UH we 
derive the E2 and Ml hyperfme selection rules, describe 
the calculations of the E2 and Ml matrix elements, and 
predict the signal ratio for relevant hyperfme transitions. 
The experimental setup is described in Sec. IIII1 and the 
experimental results and analysis in Sec. IIVI Section [V] 
comprises concluding remarks. 



II. THEORY 



(1), M = i(L + .gsS) is the non-relativistic magnetic 
dipole operator, where L and S are the spin and orbital 
angular momentum, respectively, and g s — 2.00232 is 
the electron g-factor. Notice that the expressions for the 
magnetic moments differ in the Gaussian and the SI sys- 
tems of electromagnetic units by a factor of the speed of 
light. In the conversion to atomic units this translates 
into a factor of a — 1/137. Below we use the SI system 
of electromagnetic units, in which the Bohr magneton 
is fiB = efi/(2m e ), or 1/2 atomic unit. The operator 
Q = er 2 C^ in Eq.(2) is the electric quadrupole opera- 
tor, where is a spherical harmonic tensor and r is 
the radial coordinate of the electron. 



The Hamiltonians for magnetic dipole and electric 
quadrupole transitions may be written in atomic units 
(a.u.) as (HQl 



H M1 = M B, and 
H E2 = E-Q-k, 



(1) 
(2) 



where E and B are the electric and magnetic field am- 
plitudes, respectively, and k is the wavevector. In Eq. 



We focus here on the situation relevant to our experi- 
ment, transitions betweens specific hyperfme levels F and 
magnetic sublevels to: 5p 3 / 2 (F' ,m') — ► 8pi/ 2 (F", to"). 
The hyperfme selection rules for Ml transitions can be 
derived by applying the Wigner-Eckart theorem to the 
corresponding matrix elements for a spherical component 
q of the magnetic dipole operator M, then uncoupling 
the nuclear spin angular momentum from the electronic 
angular momentum: 



J 



(8 Pl/2 F"m"\M q \5p 3/2 F'm') = a(^ il , J (2F' + 1)(2F" + 1) j^,' J,', <8p 1/2 ||M||5p 3/2 ) 



(3) 



r 



Here a is a phase factor, a = (— 1) F ~ m , the term in 
parentheses is a Wigner 3-j symbol, and the term in curly 
brackets is a 6-j symbol. The selection rules 



AF 
Ato 



and the triangle rule 



0,±1, 
0,±1, 



F' + F" > 1 



(4) 
(5) 



(6) 



follow from the 3-j symbol and the vector nature of M, 
which constrains q to the values q — 0, ±1. Repeating 
this procedure for the case of E2 transitions yields a sim- 
ilar expression for the matrix element of the g'th compo- 



nent of the spherical harmonic tensor C q 
quadrupole operator Q, 



(2) 



in the electric 



(8pi/ 2 F"m"\ Q q \5p 3 / 2 F'm') = a 



F" 2 F' 
-to" q to' 



(2F' + 1)(2F" + 1) 



J' F' I 
F" J" 2 



8pi/ 2 ||Q||5p 3 /: 



(7) 



r 



This result yields the E2 selection rules 



AF 
Am 



and the triangle rule 



0,±1,±2, 
0,±1,±2, 



F' + F" > 2. 



If the electric field and the initial atomic state are polar- 
ized, there are additional constraints on Ato due to the 

(8) double dot- product in Eq. (2). This can also constrain 

(9) the allowed values of F" under some circumstances. 

Thus, starting in an F'—O state, an Ml transition to 
F"—2 is forbidden because of Eq. (4), while an E2 tran- 
(10) sition to F"=l is forbidden because of Eq. (10). We 
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FIG. 1: Hyperfine transitions used to separate the Ml and 
E2 contributions to the 87 Rb 5p 3/2 {F') -> 8p 1/2 (F") transi- 
tion at 587.6 nm. Starting in F'=Q, the transitions to F"=l 
(dashed line) and F"=2 (solid line) are Ml allowed and E2 
allowed, respectively. The indicated hyperfine splittings are 
not drawn to scale 



TABLE I: Allowed 5p 3/2 {F') -> 8p 1/2 (F") transitions in 
87 Rb. Note that our use of F'=0 allows complete separation 
of Ml and E2 processes as shown in the first row (boldface). 





F"=l 


F"=2 


F'=0 


Ml 


E2 


F'=l 


M1+E2 


M1+E2 


F'=2 


M1+E2 


M1+E2 


F'=3 


E2 


M1+E2 



thereby obtain the key result, shown in Fig. 1, that 
5p 3 / 2 (T'=0) -> 8pi/ 2 (F"=l) is solely an Ml transi- 
tion, while 5p 3/2 (F'=0) -> 8p 1/2 {F"=2) is only E2 al- 
lowed. The various allowed transitions for 5p 3 / 2 (F') —* 
8pin(F' r ) are shown more generally in Table I for an 
unpolarizcd sample. 

In the experiment, we measure signals for specific 
F' — > F" transitions. If the initial 5p state is prepared 
without orientation or alignment, which should be a good 
approximation under our conditions as we discuss in Sec- 
tion IV, the transition rates for Ml and E2 processes will 
be proportional to the Einstein B coefficients. They can 
be obtained (l5j by summing over m" for a given mf and 
then averaging over m': 



Bmi = 



4tt2 v 
3c 2 f^, 

m ,m 



\(8p 1/2 F"m"\ M \5p 3/2 F'm')\ 
2F' + 1 



(11) 



2 1.2 



B 



E2 



TT 2 k 



E 



\5p 3/2 F'm') 



2F' 



1 



(12) 

For the special case of an F'—O initial state, used for most 
of the results reported here, the averaging is redundant as 
there is only one magnetic sublevel and the state cannot 
support either orientation or alignment. In this case the 



Ml (E2) transition strength is associated only with the 
transition to F"=l (F"=2) and the above expressions 
can be simplified to 



B 



Ml 



-BE2 = 



7T 
TT 2 k 2 

300 



>l/ 2 ||M||5p 3 /2) 

8pi/2||Q||5p 3 /2 



(13) 
(14) 



The ratio R — Bmi/Be2 is measured in the experiment 
by determining the ratio of signal sizes at a fixed laser 
intensity. Assuming that the E2 electric quadrupole mo- 
ment matrix element is known, we can determine the Ml 
magnetic dipole moment matrix element. 

Non-relativistically, the Ml matrix element between 
the 5p and 8p states vanishes identically, because the 
non-relativistic Ml operator does not couple the radial 
motion of the electrons and the non-relativistic radial 
wavefunctions of the bp and 8p states are orthogonal 
(An=0 selection rule). In other words, the entire value 
of the Ml 5^3/2 — > 8pi/2 matrix element is due to rel- 
ativistic effects. Qualitatively, this may be understood 
from the fact that the radial wavefunctions for the P\i<x 
and p 3 /2 angular symmetries arise from solving separate 
radial Dirac equations and the resulting 5p3/ 2 and 8pi/ 2 
radial orbitals are no longer orthogonal. A general rela- 
tivistic treatment of multipolar transitions may be found, 
for example, in Ref. [15]. Unfortunately, there are errors 
in the relevant formulae of that book and here we present 
revised expressions for the single-particle reduced matrix 
elements for the Ml and E2 operators: 



(nV'||M||nV 



2a 



(-k"\\C^\\k') x 



\Gn" k" (r) F n , K , (r) + F n „ K „ (r) G n , K , (r)) rdr, (15) 

(n'V'||Q||nV) = {n"\\C i2) \\K') x 
(Gn"K" (r) G n , Kl (r) + F n " K " (r) F n , K , (r)) r 2 dr.{\&) 



Here are the normalized spherical harmonics, G{F) 
are the large (small) radial components of the Dirac bi- 
spinor, and the angular quantum number k is defined 
as k = (I — j)(2j + 1). Both expressions are given in the 
long-wavelength approximation (transition wavelength is 
much larger than the atomic size). In addition, the E2 
operator is given in the conventional length gauge, which 
is more stable in calculations and does not require in- 
troducing so-called derivative terms into the many-body 
calculations. 

We evaluate the matrix elements using two meth- 
ods: the Dirac-Hartree-Fock (DHF) approximation and a 
more elaborate Singles-Doubles (SD) method [ijj ■ Both 
calculations are ab initio relativistic. The SD method 
builds upon an expansion of the atomic many-body wave- 
function into a large series of the Slater determinants 
classified by the number of excited electrons (two at most 
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TABLE II: Values for Ml and E2 reduced matrix elements, 
in atomic units. 



(a) 



Rb + 



(b) 



Rb + 



|(5p 3/2 ||M|[8p 1/2 )| 


(5P3/2||Q||8pi/ 2 ) 


DHF 


1.95 x 1(T 3 


7.05 


SD 


2.16 x 1(T 3 


6.53 



for the SD) from the reference DHF state. Major corre- 
lation corrections (e.g., core polarization) are included. 
Previously, the SD method has been shown [l7| to yield 
theoretical El matrix elements for principal transitions 
in Rb with a sub-1% accuracy. While more sophisticated 
methods are available this accuracy is sufficient for 
the goals of this paper. Numerical evaluation has been 
carried out using the dual-kinetic-balance basis set fli| 
generated in a large, 100 a.u., cavity. SD runs employed 
66 out of 75 basis functions for each partial wave up to 
£ = 5. We find that the DHF removal energies are off 
by 5% for the 5p 3 / 2 and 2% for the 8pi/2 states. Fur- 
ther application of the SD method (i.e., accounting for 
correlations) brings the theoretical energies into a 0.5% 
and 0.1% agreement with experimental values. Finally, 
the Ml and E2 reduced matrix elements are listed in Ta- 
ble II. We find that these matrix elements are largely 
insensitive to correlations (the Ml (E2) matrix element 
is affected at the 11% (8%) level). Considering that the 
SD method accounts for the dominant many-body ef- 
fects, this indicates that the theoretical accuracy of the 
SD matrix elements is at the level of 1% or better. 



III. EXPERIMENT 

The experiment is performed in a diode-laser-based 
vapor-cell magneto-optical trap (MOT), which is used 
to trap ^10 7 87 Rb atoms at a temperature of ~100 /xK 
with a density up to 10 11 cm -3 [l[. As in a stan- 
dard MOT, the trapping light is tuned 13 MHz below 
the 5si/ 2 (-F=2) — > 5p 3 / 2 (F'=3) cycling transition (see 
Fig. 2(a)). However, the repumping light, which pre- 
vents accumulation of population in 5s;l/ 2 (F=1), is tuned 
midway between the 5s 1 / 2 (-F 1 =l) - * 5p3/ 2 (-F'=l) an d 
5si/2(F=l) — > 5p3/ 2 (F'=2) transitions. This allows us 
to lock this laser to the corresponding crossover reso- 
nance in the saturated absorption spectrum. More im- 
portantly, it enables the subsequent population trans- 
fer described below to be done conveniently with readily 
available acousto-optical modulators (AOMs). This sub- 
optimal MOT configuration results in only ^20% fewer 
atoms than the standard tuning of the repumping laser 
to exact resonance with 5si/2(-F=l) — > 5p 3 / 2 (F'=l or 2). 

To spectroscopically separate the Ml and E2 contri- 
butions to the 5p3/ 2 — > 8pi/2 transition, the atoms must 
be prepared in the F'=0 state of 5p 3 / 2 - This is done 
with the excitation steps shown in Fig. 2, with the 
timing sequence shown in Fig. 3. First, the trapping 
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FIG. 2: Level diagrams showing the transitions used for state 
preparation and 5p 3 / 2 (i ? ') — 8pi/ 2 (.F") spectroscopy, (a) The 
trap (solid) and repump (dashed-dotted) lasers provide an 
ultracold sample of 87 Rb atoms which are optically pumped 
into 5si/ 2 (i ? =l) using the depletion laser (dashed) and spon- 
taneous emission (dotted), (b) The preparation laser (either 
pulsed or cw, dashed) provides excitation to 5p3/2(F'=Q), and 
the tunable dye laser at 587.6 nm (solid) further excites to 
8pi/ 2 (.F"=l,2). Atoms in 8pi/ 2 are detected by pulsed pho- 
toionization at 1064 nm (dotted). 



and repumping beams are switched off with AOMs, leav- 
ing the atoms in 5si / 2 (-F=2). Next, a depletion beam 
tuned to bsi/2(F—2) — ► 5p3/ 2 (.F'=2) is switched on in or- 
der to optically pump the atoms into 5s 1 / 2 (-F=l). This 
beam is derived from the trapping laser by frequency 
shifting with an AOM. Finally, a preparation beam, de- 
rived from the repumping laser, is switched on to ex- 
cite the 5si/ 2 (F=l) — > 5p3/ 2 (-F'=Q) transition. The 
5^3/2 8.P1/2 transition is excited with 587.6 nm light 
from a cw single-frequency tunable ring dye laser system 
(Coherent 699-29 with Rhodamine 6G dye) pumped by 
an argon-ion laser. This light is left on throughout the 
excitation cycle. 

Although the 5s!/ 2 (F=l) — » 5p 3 / 2 (F'=0) transition 
is a cycling transition with respect to F (i.e., F'=0 can 
only decay back to F—l), the F'=0 population is ex- 
pected to be transient due to optical pumping. Be- 
cause the 5si/ 2 (F=l) level has a higher degeneracy than 
the 5p 3 / 2 (F'=0) level, after a few cycles of excitation 
and spontaneous decay, most atoms are pumped into 
dark magnetic sub-levels and cannot be further excited. 
If the excitation light is 7r-polarized and drives Am—0 
transitions, the population accumulates in 5s 1 / 2 (F=l, 
m=±l). In order to determine the optimum timing for 
the various beams in the presence of this optical pump- 
ing, we solve the rate equations for the time-dependent 
level populations of the system. The 5s!/ 2 (F=l,2), 
5P3/ 2 (^'=0,1,2,3), and 8p 1/2 (F"=2) levels are included. 
In addition, 5si/ 2 (F=l, m=±l) is included as a dark 
state that can be populated by spontaneous decay, but 
not excited. Appropriate lifetimes, transition strengths 
and branching ratios are used in the simulation, as well as 
the experimental values for intensities and timings. The 
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FIG. 3: Timing diagram for the 5p 3/2 (F') - 8p 1/2 {F") spec- 
troscopy. The case of pulsed preparation (5si/2 — 5^3/2 exci- 
tation) is shown. Populations in the 5p3/2(-F'=0) (solid line) 
and 8pi/2(F"=2) (dashed line) levels, calculated by rate equa- 
tions, are also shown. The preparation and excitation inten- 
sities are 5.2 mW/cm 2 and 62 W/cm 2 , respectively. 



long-lived 8pi/ 2 level (480 ns lifetime [III) is assumed to 
decay back to 5p 3 / 2 via 6si/ 2 . The resulting 5p 3 / 2 (F'=0) 
and 8pi/ 2 {F'=2) populations, shown in Fig. 3, reveal an 
optimum preparation time before optical pumping into 
the dark state occurs. However, as discussed in Sec. IIVI 
this dark-state pumping is not as severe in the experi- 
ment as these calculations predict. 

For detection of 8px/2 atoms, we use photoioniza- 
tion with an injection-seeded pulsed Nd:YAG at 1064 
nm. The fluence of each 5 ns FWHM pulse is typically 
0.11 J/cm 2 , yielding an estimated 8pin ionization prob- 
ability of 63%. The entire excitation sequence is syn- 
chronized with the 10 Hz repetition rate of this pulsed 
detection laser. We note that the 532 nm second har- 
monic on the Nd:YAG laser was originally used for 8pi/ 2 
ionization, but this yielded a significant background due 
to two-photon 5si/2 ionization. Using the 1064 nm light 
mitigates this problem, since 5si/ 2 ionization now re- 
quires four photons. Also, lower pulse energies can be 
used because the 1064 nm photoionization cross section, 
1.62 xlO -18 cm 2 (calculated using phase-shifted Coulomb 
wave functions), is a factor of 6.4 larger than that at 532 
nm. 

The MOT is located between a pair of 95% transpar- 
ent grids separated by 2.09 cm. A pulsed electric field of 
~216 V/cm is applied about 1.1 /is after the photoion- 
ization pulse in order to extract the resulting photoions. 
These ions are detected with a discrete dynode electron 
multiplier (ETP model 14150) whose output is sent to 
a boxcar averager in order to select the desired time-of- 
fllght window. A computer records the data from the 



averager as the 5p 3 / 2 — > 8pi/ 2 laser is scanned. 



IV. RESULTS 

A typical spectrum for the 5p 3 / 2 (F'=0) — > 8pi/ 2 
(F"=l,2) transitions is shown in Fig. 4(a). The 
5p 3 / 2 (F'=0) -> 8p 1 / 2 (F"=2) E2-allowed transition is 
clearly seen, while the 5p 3/2 (F'=0) -> 8p 1/2 {F"=l) Ml- 
allowed transition is not observable above the noise. This 
allows us to set a limit on the Ml matrix element. 
For comparison, we show in Fig. 4(b) a scan over the 
5pa/2(F'=2) -» 8p 1/2 (F"=l,2) transitions, both of which 
are E2 allowed. Since the depletion light transiently pop- 
ulates 5p 3 / 2 (F'=2), and the 8pi/ 2 excitation light is al- 
ways on, detected ions are associated with transitions 
from this level. This scan is useful on two accounts. 
First, it allows us to check the expected ratio of m- 
averaged transition strengths for these two E2-allowcd 
transitions. We measure a ratio of areas under these 
peaks of 0.43±0.04, which compares well with the pre- 
dicted value of 0.43. This prediction is based on a ho- 
mogeneous distribution of population in the magnetic 
sublevels, a reasonable assumption in the complex ra- 
diative environment of the MOT. Second, it provides 
an accurate calibration of the 64.24 MHz 8p 1/2 (F"=l,2) 
splitting [2(J, thereby telling us where in Fig. 4(a) 
the 5p 3/2 (F'=0) -» 8p 1/2 (F"=l) Ml-allowed transition 
should be located. 

Due to nonlinearity in the scanning of the dye laser, an 
uncertainty in this frequency calibration on the order of a 
few percent is expected. We verify that this uncertainty 
is <5% by comparing scans over the 5p 3 / 2 (F'—2) — > 
8p 1/2 (F"=l,2) and 5p 3/2 (F'=3) -f 8p 1/2 (F"=l,2) tran- 
sitions. The latter is obtained by leaving the MOT beams 
on and keeping the depletion and preparation beams off, 
which results in populating the 5p 3 / 2 (F'—3) level. The 
uncertainty in the 8pi/ 2 (F"=l,2) splitting due to scan 
nonlinearity is taken into account in the fits that are used 
to set a limit on the Ml transition. 

As discussed in Sec. IIII1 we expect that optical pump- 
ing into the 5si/ 2 (-F=l) dark state will result in an op- 
timal timing of the 5si/ 2 (F=l) — > 5p 3 / 2 (F'=0) prepara- 
tion pulse. However, we find that a cw preparation beam 
actually yields a larger signal (by a factor of ~2) and a 
better signal-to-noise than a pulsed preparation beam. 
We believe that this counter-intuitive suppression of the 
optical pumping is likely caused by multiple scattering of 
resonance fluorescence light in the optically dense sample 
and/or Zeeman precession of the spin in the quadrupole 
field of the MOT 

To set a limit on the strength of the Ml transition, 
we fit the observed 5p 3/2 (F'=0) -> 8p 1/2 (F"=2) E2 sig- 
nal and the unobserved 5p 3/2 {F'=0) -> 8pi/ 2 (F"=2) Ml 
signal to a pair of Gaussians separated by the 64 MHz 
8p 1 / 2 (_F"=l,2) splitting. Although we do not understand 
the lineshape of the observed transition in detail, the 
scans fit quite well to Gaussians with a typical full width 
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FIG. 4: (a) Solid line: Typical spectrum of the 5p 3 / 2 (F'=0) 
— > 8pi/ 2 (F" ~1,2) transition. As shown in the inset, no signal 
is observed for F"=l (near -64 MHz), setting a limit on the 
relative Ml transition strength. Dashed line: Best fit to a 
pair of Gaussians, as described in text. Dotted line: residuals 
from the fit, offset by seven units for visibility, (b) Spectrum 
of the 5p 3/2 (-F'=2) -> 8p 1/2 (F"=l,2) transition (solid line), 
with best fit to a pair of Gaussians (dashed line), and residuals 
offset by seven units (dotted line). In this case, allowed E2 
transitions to both F"—l and 2 are observed. 



at half-maximum (FWHM) of 10 MHz, as seen in Fig. 
4(a). The 480 ns [19| lifetime of the 8p±/2 upper state 
of this two-photon transition would yield a Lorentzian 
natural linewidth of 330 kHz. However, the two lasers 
driving this transition each have linewidths in the range 
of 1-2 MHz. Also, the inhomogeneous magnetic field sam- 
pled by the atoms in the MOT results in a distribution 
of Zeeman shifts for the various magnetic sublevels. The 
fact that the first step of the two-photon transition is on 
resonance and close to saturation may also affect the line- 
shape. The second step (E2) transition is very weak, with 
a transition rate estimated to be ~10 4 times smaller than 
that of the first step (El) transition. The fits are clearly 
not consistent with a Lorentzian lineshape with the 5.9 
MHz FWHM natural linewidth of the 5pz/2 intermedi- 
ate state. Since the observed peaks for other transitions, 
such as those shown in Fig. 4(b), are also Gaussian with 
similar widths, we feel justified in using a Gaussian fit for 
setting a limit on the unobserved Ml transition. Under- 
standing the details of the lineshape and its dependence 
on various parameters is an interesting topic for further 



investigation. 

An upper limit on the Ml transition strength is deter- 
mined by least-squares fitting of a number of individual 
scans to a pair of Gaussians. The separation between the 
peaks is constrained to 64 ± 3 MHz, and the FWHM of 
each Gaussian is allowed to vary independently from 7 
to 21 MHz. We also include constant and linear terms 
in the fits to account for the background and its slow 
drift, respectively. A total of 32 scans are included, 20 of 
them using pulsed preparation and the remainder using 
cw preparation. An example of the latter is shown in 
Fig. 4(a). For each scan we determine the ratio R of am- 
plitudes for the two Gaussians, for comparison with the 
theoretical prediction obtained using Eqs. (13-14) and 
the matrix elements in Table II. The mean experimental 
ratio is found to be i?=0.0026 and the standard deviation 
of the mean is 0.0020. Assuming a normal distribution 
of values, we set an upper limit on R by adding to the 
mean two standard deviations of the mean. This yields 
R <0.007 with 95% confidence. By this criterion, our 
result is consistent with zero. We verify that omitting 
the linear background term from the fits has a minimal 
effect, changing the mean value of R to 0.0024 and its 
upper limit to 0.006. We perform a similar statistical 
analysis of the scans by simply integrating the signal in 
6.4 MHz wide bins instead of doing Gaussian fits. This 
yields a mean value of i?=0.0021 and an upper limit of 
0.005. As our final value, we quote R <0.007. 

Statistics connected to an F-test [22| reveal that the 
Gaussians fitted to the position of the Ml peak are not 
significant, i.e., the probability of an actual peak at the 
expected location is well below 95%. When we add a 
simulated peak to our data, the F-test shows that, with 
our present noise level, we need to add a signal with 
R >0.01 in order for it to be significant. 

In Table III we summarize our results. In our experi- 
ment, the limiting value of the ratio R of signals for the 
5p3/2(^'=0) -> 8pi/ 2 (F"=2) Ml-allowed transition and 
the 5p 3/2 (F'=0) -> 8p 1/2 (F"=2) E2-allowed transition 
is determined directly. The previous polarization-based 
measurements Q can be converted to yield an indirect 
value for i?, shown in the middle column of Table III. 
Obviously, our upper limit on R is well below the range 
reported in this earlier work. Also shown is our theoreti- 
cal prediction. The experimental limit is consistent with 
this theoretical prediction, but we are not able to pro- 
vide a statistically significant non-zero result with which 
to compare the theory. Table III also compares the Ml 
matrix elements resulting from our measurement with 
the earlier experiment and the theory. 



V. CONCLUSIONS 

We use high-resolution laser spectroscopy of ultracold 
87 Rb atoms to set an upper limit on the Ml contribution 
to the bp — 8p transition. Starting in an F' =0 state, we 
take advantage of the hyperfine selection rules, which al- 
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TABLE III: Signal ratios R and Ml matrix elements for 
the 5p 3/2 (F'=0) -> 8p 1/2 (F"=l) and 5p 3/2 (F'=0) -» 8p 1/2 
(F"=2) transitions. 





Based on 


Present 


Present 




Ref. [7] 


Experiment 


Theory 


M1/E2 ratio, R 


0.38 ± 0.03 


< 0.007 


0.0006 


Ml matrix 








element (10 -3 a.u.) 


55 ± 1.4 


< 7.4 


2.16 



low for complete separation of the E2 and Ml processes. 
Our upper limit on the Ml transition strength is consis- 
tent with the theoretical prediction presented here, but 
disagrees with a previous measurement based on the lin- 
ear polarization degree of the transition Q • 

Our present limit is based on the signal for the Ml- 
allowed hyperfme transition being <0.007 of that cor- 
responding to the observed E2-allowed hyperfme transi- 
tion. With an improved signal-to-noise ratio, the Ml 
transition should actually be observable, which would 
allow the theoretically predicted ratio of i?=0.0006 to 
be tested. However, we are exciting the atoms to 
&P1/2 by a two-photon process where the first step 
is resonant with the 5si/2(.F=l)— > 5p3/ 2 (F'=0) tran- 
sition. When the second step is tuned to the Ml- 
allowed 5p 3 / 2 (F'=Q) — > &P\/ 2 {F" =1) transition, the 
5s 1/2 (F=l) -» 5p 3 /2(-F"=l) -» 8p 1/2 (F"=l) two-photon 
transition, where the second step is E2-allowed, is also 
two-photon resonant, but with the F'=l intermediate 
state off resonance by 72 MHz. Accounting for this de- 
tuning and the relative hyperfme transition strengths, we 
estimate that this off-resonant two-photon E2 contribu- 



tion to the signal at the Ml resonance would give an 
apparent i?=0.0007, which is coincidentally quite close 
to the predicted Ml contribution. These contributions 
could be distinguished by varying the detuning of the 
first step of the excitation. We note that Stark mix- 
ing of s or d character into the 8p level by an external 
field could give rise to a weak El contribution at the 
location of the Ml transition. However, this mixing is 
negligible under our conditions. An alternative method 
for measuring the Ml contribution would take advantage 
of the fine-structure selection rules for Ml and E2. The 
5pi/2 ~ * 8??i/2 transition at 579.5 nm is Ml allowed but 
E2 forbidden, so a similar experiment to that presented 
here, but without the requirement of hyperfme resolu- 
tion, could be performed. Light at 795.0 nm would be 
needed to populate 5pi/2- 

Our technique for measuring Ml matrix elements can 
be extended to other states in 87 Rb and to other atomic 
systems possessing an F'=Q state. For example, the ra- 
dioactive isotopes 223 Fr and 225 Fr both have nuclear spin 
1=3/2 and therefore an F'=0 level in their 7p 3 / 2 state. 
Since relativistic effects increase with atomic number Z, 
heavier atoms such as francium will have larger Ml ma- 
trix elements, thus serving as good test cases for atomic 
structure calculations. 
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